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We consider B — > PS decays where P stands for pseudoscalar and S for a heavy ( 1500 MeV) scalar 
meson. We achieve agreement with available experimental data - which includes a two orders of 
magnitude hierarchy - assuming the scalars mesons are two quark states. The contribution of the 
dipolar penguin operator On is quantified. 

l ■ I. INTRODUCTION 

The scalar sector below two GeV is poorly understood, nevertheless several features -like the presence of two 
multiplets and several of their properties - naturally arise in the analysis of a number of authors. A first set 
of scalars with masses around 1.5 GeV [l[ are grouped in a heavy multiplct, including the Kq (1430), ao(1450), 
/o(1500) for the octet, /o(1370) which is identified with the singlet and the /o(1710) which seems to be mainly 
glucball. The octet is nearly degenerate, like similar pseudoscalar, vector, axial vector and tensor multiplets, 
{-^ | their widths are small (< 100 MeV). The mixing angles seems to be small except by the singlet-glueball which is 
(f) . around —20°, according to H. Y. Cheng in ref. 0j. It has been more difficult to establish the lighter multiplet, 
C\| ' even the existence and nature of some of their members is in doubt. The light multiplet should include the 
qq , ao(980), /o(980) and the k = A'q(800) in the octet; while the singlet could be identified with the a = /o(600). 

The mixing is not clear and their widths are very large. Ideally, the former multiplet can be identified as the 
ground state of quark antiquark bound states with angular momenta one while the later with the ground state 
of four quarks systems with zero angular momenta. In the real world an undetermined mixing between the two 
multiplets is expected. Alternatively both multiplets could be identified as quark-antiquark states with angular 
momenta one, the lighter being the ground state while the heavier the first excited state. 
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The full understanding of the scalar multiplets previously described remain a challenge, both from the exper- 
imental perspective as well as from the theoretical point of view [l|. To start with, there is not enough and 
conclusive experimental information regarding the existence and properties of the scalars. Notice that the infor- 
mation is poor not because of the lack of sources of scalar mesons, for example many of the decays of particles 
containing c or b quarks involves the production of scalar mesons. The information on the scalars is scarce 
because of the large width they have since that produces a large overlap with nearby resonances and with the 
background. In spite of those problems, precise experimental results are available |l|, [H, 0| for the mass and 
width of the /o and Kq, for the (3 angle Q of the CKM matrix and for several partial widths. It has been 
speculated that the a angle can be extracted in processes involving scalars [|| and new projects like the LHCB 
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will improve the old measurements and obtain new results. Relevant to our work are the branching ratios 
for the B — > PS decays measured by different groups, which show a non trivial hierarchy. The experimental 
data collected in Table I suggest that, for B — > PS decays including members of the heavy scalar multiplct, the 
order of magnitude of the branching ratios involving the Kq (1430), the / (1370), / (1500) and the a (1450) are 
different. 

On the theoretical side the situation is not better. The origin of the difficulties are the non perturbative regime 
of QCD and the limited computer capacity for the lattice approach. The nature of the observed scalars has 
been discussed at length and proposals exist to identify them as 2 or 4 quark states, glueballs, molecules, 
etc. and several theoretical formalisms have been developed to calculate non leptonic decays. The simplest 
one is the so called 'Naive Factorization Approach' (NFA) @, which in general produces the correct order of 
magnitude and its predictions are in rough agreement with the experimental results. Discrepancies are known 
to occur in two cases, for 'color suppressed' processes and when important re-scattering effects arc involved, for 
example processes where direct CP violation is relevant 0, [l(| . The advantage of formalism where a systematic 
expansion is implemented and where higher order correction can be organized and controlled are of great 
importance (QCDF, SCET, pQCD, LCSR, etc. [UGlGSEi), 

in particular when high accuracy predictions 

are required. 

Additional reasons to study the B — > PS decays are: they offer a window to study the spectroscopy and the 
dynamics of the scalar sector, the B — > 3P decays get a contribution from the B — > PS, PV, PT, so that 
in order to achieve an appropriated estimate for the former decay the latter must be well known [15|]. In a 
similar way one can argue that in order to extract signals of possible new physics, the contribution of low lying 
conventional physics has to be known in detail, including the contributions of the scalar mesons [l6| . We believe 
that the understanding of the physical origin of the hierarchy of scales appearing in the B — > PS decays can 
shed some light on the nature of the scala rs fl7l . [iH ] . Complementary information on the nature of the scalars 
may be obtained from D — > PS physics [19|]: in the first case through the decay constants, f s while in the 
latter through the F DS form factors. The purpose of the present work is to consider the B — * PS decays 
with S a member of the heavy scalar multiplet. We assume that the leading contribution to these processes 
is given by the NFA and that, in first approximation, contributions other than the leading one can be safely 
neglected. In these conditions the dominant contribution can be clearly identified and the existence of the scales 
in the branching ratios naturally arises. Besides the NFA our approach can be summarized along the following 
lines: we include ten dimension six four quark operators and the dimension five chromomagnctic operator On 
[20| . annihilation contributions are included and the form factors required are obtained by using sum rules, 
so infrared divergences are absent. This approach, together with SU(3) symmetry, allows us to reproduce the 
pattern observed experimentally. 



II. BRANCHING RATIOS AND MIXING 



Our results are summarized in Table I. It is worth remarking that both the experimental data and our results 
points to the existence of branching ratios that ranges from 45 to 0.5 (in units of 10~ 6 ). In the following 
paragraphs we introduce the notation, conventions and explain the procedure we follow to obtain these branching 
ratios. Within the NFA the hadronic matrix elements can be reduced to products of decay constants and form 
factors. In order to achieve this one uses the 'vacuum saturation' approximation and neglect other intermediate 
states. This seems to be a reasonable assumption since the hadronic resonances have masses in the 1 — 2 GeV 



range, far from the nib region. For the invariant amplitude we write A4 



while the branching ratios are given by B = tb G % \ A \ 2 p / 1 6 irm \ 
The decay constants and form factors are defined as % O El: 
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TABLE I: Branching ratios for the B — > PS decays (in units of 10 6 ), for the heavier scalar multiplet. The values 
reported for the widths marked with * include the corresponding branching of the scalar decaying channel. To obtain 
the NFA predictions we used B(/o(1370) 2tt) = 0.26(1), B(/ (1500) -> 2tt) = 0.35(2) and for the a (1450) -> ttt? no 
reliable value exists 0|. 



with q=p 1 - p 2 . 

We have left to the appendix details regarding the effective Hamiltonian we use - which includes ten dimension 
six operators and the so called On operator - and the matrix elements evaluation. The most interesting decays 
are those involving the S = ATq(1430) both because they have the largest branching ratio (around 40, in units 
of 10~ 6 ) and because the theoretical predictions are the cleanest. The a$ term is by far the dominant one. The 
amplitudes are proportional to \t s fK*o-6 rr iK<-/ m s'>Tib ~ Xts^rtibrriK* j k- times SU(3) factors. The origin of the 
enhancement is a combination of large CKM matrix elements, the novanishing decay constant and a large ttik * 
(Chiral enhancement) mass. The SU (3) symmetry allow us to relate different decays involving the Kq and so, 
by measuring one of them, one can predict the others, a fact that is not distorted by the On contributions. For 
the numerical analysis we used the following input parameters: F B7T = 0.27(4), F BK = 0.33(4), m s (2.1) = 90 
MeV, ,F BQ o(i450) = F BK*(iAm) = Q 2 q arK j 7 when required, SU(3) relations are invoked. Although predictions 
for fx* are available [17| |. we prefered to include the B + — > 7r + A^*° experimental value as an input, obtaining 
thus ffc» ~ 58 MeV {fx* — 56 MeV when the On is taken into account). The branching ratios we obtain for 
other channels involving the Kq are reported in Table I. Notice that the value obtained for Jk* is not far from 
the theoretical predictions (see Table II). 

We now consider the decays involving S = /p(1370), /o(1500) and /o(1700). Their relevance stem from the 
large branching ratios predicted for them - of the same order as the Kq - and also due to the possible 
glucball nature of the /o(1700). Their amplitudes are proportional to XtsdembrriK* ff times SU(3) factors and 
mixing angles (s- quark content). Our predictions for these processes are included in Table I, unfortunately the 
experimental results are still inconclusive. Note that except the /o(1500) decay channel, the NFA plus SU (3) 
symmetry for the heavy scalar multiplet leads predictions for the branching ratios in rough agreement with 
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the experimental values. However, even if the experimental data is poor the discrepancy between our results 
and experimental data is evident, there is a one order of magnitude difference. In this sense it is important to 
remark that in order to obtain the results of tabler I we assumed, following H. Y. Cheng a mixing between 
the glueball, singlet and octet components given by: 
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where = sin0j and so on. The angles 812 — 2°, #13 ~ —21° and #23 — 2° are the mixing between singlet- 
octet, singlet-glueball and octet-glueball, respectively. The singlet and the octet are /o(1370) ~ /sing. = 
yplfl N + S/y/3, / (1500) ~ / oct . = N/VS - S^/2/S, S = ss, N = {uu + dd) /y/2 and G = gg the glueball. 
Thus, in this approach 0], there is only an small mixing between the singlet and the glueball. Using these values 
the prediction for B — ► /o(1500)if is in conflict with the experimental data. One way to avoid this problem is 
to leave #12 as a free parameter, keeping the others fixed. Using the experimental data we obtain the following 
inequality for the mixing between the singlet and the octet: 



Wi+c 12 j| < 0.34. (3) 



These constraints lead two possible values: 



35° < 0i2 < 74° (4) 
215° < 6<i2 < 254° (5) 

It is worth noticing that these values for the mixing are close to those mentioned by several groups [l[. 

Finally for the decays involving the ao(1450), B — > ao(1450)7r, ao(1450)if, the terms proportional to 04 — a@ ~ 
almost vanish and the branching ratios are smaller. Two different cases must be considered. The first when the 
amplitude is dominated by the tree level contribution ai (The amplitudes are proportional to \ u d a i m2 Bfir), then 
the theoretical prediction is reliable and the branchings are predicted to be are around 10 (in units of 10~ 6 ). 
The second case arises when no tree level contribution exist and terms like annihilation are dominant. In this 
case the branchings arc of order 0.1-1 (in units of 10 -6 ) but the theoretical uncertainties arc larger since other 
contributions (FSI for example (Iol|) maybe important. Unfortunately little is known about these corrections. 



III. SUMMARY 



In this work we studied the B — > PS decay where S stands for a member of the heavy scalar multiplet. The 
computation have been done assuming the heavy scalar multiplet is a two quark states, using SU(3) symmetry 
and the naive factorization approach. Our conclusions can be summarized as follows: 

• Within the error bars, it is possible to reproduce the hierarchy of branching ratios experimentally observed 
in the B — ► PS decays, whether or not the operator On is included. 
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Ref. 


(///)x*(1430), / 


/a (1450) 


// (1500) 


m s [GeV] 


Meurice-87 [21] 


27 


- 


- 




Narison-89 [21] 


40(6) 


- 


- 




Maltman [21] 


42(2) 


390(159) 


- 




Chernyak-01 [21] 


70(10) 


- 


- 




Shakin-01 [21] 


30 


207 






Pennington-01 [211 










Du-04 [211 


42(8), 427(85) 






0.14 


Cheng-05 [17J at /x = 1 GeV 


445(50) 


460(50) 


490(50) 


0.119 


Cheng-05 [17] at = 2.1 GeV 


550(60) 


570(60) 


605(60) 


0.09 


lattice-06 [1] 











TABLE II: Decay constants for scalars (in MeV). The heavy scalars are assumed to be two quark states. Notice that 
the constants computed by Cheng, were obtained by using sum rules, OPE and Renormalization Group equations that 
render / scale dependant. 

• When the singlet-octet mixing given by [l[ is used, we obtain a prediction for the /o(1500) which is 
one order of magnitude above the experimental limit. A solution to this problem can be obtained by 
modifying the mixing matrix. In such a case one obtain a constrain on the singlct-octet mixing and its 
s-quark content. 

• The contribution of the On operator is around 30 % in decay channels involving the Kq. The On 
contributions approximately keep the SU(3) relations between different decay channels. 

• The chiral enhancement predicted by the NFA could be used to test the quark structure of the heavy 
multiplet. Strong deviations from the NFA results could be interpreted as a signal that the heavy scalars 
are not pure two quark state. 
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APPENDIX: NAIVE FACTORIZATION APPROACH (NFA) 

The relevant effective Hamiltonian is given by [1] : 



We// 



G p 
V2 



10 



v ub v; q (do? + c 2 o u 2 ) - v tb v; q (J2 C *°* + c 9°9 



vi=3 



+ h.c. 



(A.l) 



where \' q = V q 'bV*,„, with q = d,s, while q' = u,c,t. The Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements are denoted by Vij . Oi stand for the following four fcrmion operators: 



Ox = (qu) L {ub)! 



02 = {u a bf3) L {qi3U a )l 



O3 = (qb)L E,'(«V)ii 04 - {q a bp) L Y,%l' a )L, 

q' 

05 = (qb) L E q '(q'<l')R, 06 = (q a bp) L J2%^)n = ^^(^S-P^S+P. 

9' 1' 

07 = |(#)LE,.e,,.(9'g')j!, 08 = ^(q a b/3)L^2e ql {q' fj q' a ) R = -3^2e q ,{q'b) s „p(qq') s+P , 

q> q> 

09 = l(qb) L E q ' e q >(<i'q')L, o w = |fe^) L ^ eg ,(^^) L = f ^e^kfe')^ 



(A.2) 



and the dipole penguin operator: 



0ii = ^m b qa^RT a bG^; 



with T a = A a /2 the SU(3)c generators. The Wilson coefficients Ci appear in the combinations ct2i— 1 = 
C 2 i-i + C21/N, a 2 i = C*2i + C21-1/N. The numerical values are taken from [8[. Similarly we define an = 
(8/9)a, s Cii (m 2 /47ra 2 ) ~ —5.7 • 10~ 3 . For the gluon momentum we use q~p b —p s ~p B — Pk/2, so g 2 ~ m 2 B /2 
[2(j |. Taking a s (g 2 ~ 7713/2) = 0.21 and Cn = —0.29 one obtains an ~ —5.7 ■ 10~ 3 . Chiral projections are 
L, R = 1 =F 75 . Using the relation 2(Ti) a p(Ti) 1 $ = <5 Q 5<5^ 7 — {1/N)8 ai 8f3s and the Fiertz reordering one obtains 

for On Hi 



.Gf. C g asm b 



" 1 r 2 a a ' 

iV 2 aP 75 " Nc Ta ^ & 



fc M [3i«a-Ry^^« - SiqaRq'pq'^Rbs + q a R lp q' p q'y^ Rb S - q a R^ v ^^Rbs] (A.3) 
where A: 2 ~ tob/2 — m 2 K /%. 

The amplitudes, including On contribution, in the NFA are given by: 



j4ro 



A„ d (aiX| 0a+ +a 2 X^ +7r _^ ) - A 



(a 6 + a 8 )m 2 • 
04 + aio — — ; 7 I A so + 



m(mb + m u ) J B a o 



+ 2(a 3 - a 5 ) + a 4 + 



a 9 - a 7 - arc (as - a 8 /2)m 2 



m u (mb + rrid) J (a o w ) 



*7 + 



8 



A 



B°- 



Kd(aiX^' ^ + + a 2 X ( ?_ w+ ^) - X td 
+ ( 2{a 3 - a 5 ) + a 4 + 



(a 4 + a 10 )Xg° On+ - 2(a 6 + a 8 )Xg° ^ 



a-9 - a 7 - a w (a 6 - a 8 /2)m| \_ ^ga 



m u (m b + m d ) ) K 77 ) 



^ud,0>l{Xg- go + X sayr -) — X t d 
+ ^a 4 + a 1Q 

And 



a 4 + aio - 



* 1 



rh(m b + vn u ) 



(0 6 + Og)m| \ V B~ , r,„^v3 U d 



m(m b + m u ) 
^B%o+X^ ) + a 2 X< 



X$> v - + (a s - 2a 6 )X»± n _ 



Hd 



(a 4 + aio) ^s-^o - 2(a 6 + a 8 )X°°_ 7r0 



3 1 (a e + a 8 )ml \ ^ / (a 6 + a 8 )m| i 

- ( a 4 - -(o g - a 7 ) - -aio - 



m u (m b + m d ) J ^ B a + \ 4 + 010 »5i(ra(, + m„) / "^ a o 7r ° 



.4 



3 (2a 6 - a 8 )mj \ w ° 

^ a 7 + 7; 7 C I A B°a° 



2m d (m b - m d ) ) B a o 



9 3 (2a 6 - a 8 )m| \ 



(a 4 + aio - (as + a 8 )r^X§ 0a+ 



«4 — 2 fll ° 



(2a 6 - a 8 )m| s x& 
(m b + m d )(m s + m d ) ' 



Ab~ —>K- S° — ^us a l 



^B-S° + ^SOif- 



(a 4 + ai ~ {af- + a 8 )rf ) X*_ s0 



+ I O4 + OlQ - 



2(a 6 + a%)m\ 



(m b + m u )(m s + m u ) 



X 



S°K- 



A 



B~^K°a' 



^us a lX 



K°a- 



+ a 4 + aio 



(04 - \a w - (ae - a 8 /2)rf )xf! a _ 



2(a 6 + a 8 )m\ 



(m b + m u )(m s + m u ) J K ° a o 



X* 



2a, 



6 — a 8 



30 

—an j X B _ 



A^ 



B a ^K»S° 



30 

~ 32 ' 5°k° 



a 4 



aio (2a 6 - a s )m 



2 (m b + m d )(m s + m d ) 



A ; 



Ab°^tt+kz- — ^usdlX = C _ + - A 



(a 4 + am - (af- + a 8 )r*) 



r* \X K ° 
'xr^B°iv+ 



9 



+ ( a 4 - -a w 



(2a 6 - a 8 )m 



(m b + m d ){m s + m d ) 



X' 



^■B-^-k-K* 



^us^lX — \ts 



a 4 - ioio - (of - - a 8 /2)r*)x*°°_ 



.4 



«1 



■4 + aio - 

+ | 4 + aio 
A 



2(a 6 + a 8 )m| \ 
(m b + m u )(m s + m u ) J R o"^' 



a 2 X 



B-K* 



A, 



(a 4 + aio - (af- + a 8 )r*) X*l~ 



XI y ^B-ir° 



2m 2 B (a 6 + a s ) \ B - 3 
(m b + mu )(m s + m u )) + ^ ~ a ^ X B~K- 



a ^ X B U -K*- 



- A* 



(o4 + oio-(af-+a8)r*)X^ 



[ |(fl9 - a 7 ) - 2a 4 + aio + (2a 6 - a 8 )r* + ^ (12 - 30r* - l) J 



+ a 4 + aio 



2m 2 B (a 6 + a s ) 



(m b + m u )(m s + m u ) J K o 71 



X, 



Xus^X^g — A 



a 4 



aio (2a 6 - a 8 )m| 



2 (m b + m d )(m s + m d ) 



a4 -^-(a§ ff --ag/2)r*)xi 



Xj B°rj 

an 



+ ( 7r( tt 9 - or) - 2 «4 + aio + (2a 6 - a 8 )r* + (12 - 30r* - 1) ) X% 



B°K*° 



aio (2a 6 - a 8 )m| \ B o 

a 4 — — n 7 w \ I -^i?* 

2 (m b + m d )(m s + m d ) J "o v 



(A.4) 



and An 



B°^K 



„M 



» +7r _ = 0, 5° = ao, c and / , r* ~ 2m 2 K ,/m b m s , r^f ~ m 2 /m b m s and tJ 1 = a 6 r* 



a n [I2(l-rf)-l] /32. 



X£_ /o =< K-\(us) L \0 X f \(ub) L \B- >= f K (m 2 B - m 2 a )Ff fo {ml) = / x '» ^"" K/ (m| - m^) shnfe 



^Ic°/ =< ^°I(^)l|0 X / |(<&) £ |B° >= /*o(m| - ro? )f * f °(m%) = ^ K 1 (m% - mjjsh^s 
X§~, =< K*-\(su) L \0 X n + \(ub) L \B Q >= f K *-(m B - ml)F B °* + (m 2 K ,) = -f K -{m\ - ml)F B (m 2 K ,) 



XbK- =< K* o »\(sd) L \0 >< ir-\(db) L \B- >= f R »o(m B - m*)Ff * «,) = -/*.«. (m£ - mi)F*~* (m 2 K >) 



X;^ a =< K*-\(su) L \0 >< n Q \(ub) L \B- >= f K -(m B - m^)F Q fl " T «K;) = -/x, 



10 



X f 0R , =< A(uu) L \0 >< K* °\(sb) L \B° >= A( ro | - ^1^^) = A (m | _ m ^ )r ^ F f Qo+ (m^) 



_^ 7^*01 Irw^„0|/Ji\ inO^_ f / 2 2 \pflV/„ 2 \_r / 2 ,2 \ ^(f 77 ( TO *i 



*f£o =< /C|( S -d) L |0 >< A(db) L \B° >= f K ,o(m% - ml )F s °*°(m 2 K ,) = f K . (m| - ro*) 



x/2 



A' 



2 ™2 „,2 ™2 



*£-V =< /ol-|0 >< W >= ^^^^"'"K) = ^3^^m A /^-K) 
*£feo =< fo\MO >< K°\sb\B° >= m /,// G ^rK) = (A.5) 
for annihilation: 

X foA- =< /°^"l(*«)£|0 >< 0\(ub) L \B- >= -f B (m 2 s - m K )Ff K ~(m 2 B ) 
X f<o fo =< K°f \(sd) L \0 X 0\(db) L \B° >= -f B (m% m 2 K )F^ R "( m 2 B ) = X%Z fo 
=< if*-7r+|(sd) L |0 >< 0\(db) L \B° >= -f B (m 2 K , - ml)Ff°~** ' (m%) 
=< KS°7r-|(Su) £ |0 X 0\(ub) L \B- >= -f B {m 2 K , - ml)F«°° v ~ (m|) 
X*:^ =< K*-7r°\(su) L \0 >< 0\(ub) L \B- >= -f B {m 2 K , - m^F^'^ (m%) 
X R^o =< KZ\°\(M) L \0 X 0\(db) L \B° >= -f B {m 2 K , - ml)F^\m 2 B ) 

with r I<7r = F BK /F B « ~ f K /U ~ 1.21(9) 



(A.6) 



